Back-action evasion and squeezing of a mechanical resonator using a cavity detector 
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We study the quantum measurement of a cantilever using a parametrically-coupled electromag- 
netic cavity which is driven at the two sidebands corresponding to the mechanical motion. This 
scheme, originally due to Braginsky et al. [V. Braginsky, Y. I. Vorontsov, and K. P. Thorne, Sci- 
ence 209, 547 (1980)], allows a back-action free measurement of one quadrature of the cantilever's 
motion, and hence the possibility of generating a squeezed state. We present a complete quantum 
theory of this system, and derive simple conditions on when the quantum limit on the added noise 
can be surpassed. We also study the conditional dynamics of the measurement, and discuss how 
such a scheme (when coupled with feedback) can be used to generate and detect squeezed states 
of the oscillator. Our results are relevant to experiments in optomechanics, and to experiments in 
quantum electromechanics employing stripline resonators coupled to mechanical resonators. 

PACS numbers: 



I. INTRODUCTION 

Considerable effort has been devoted recently to at- 
tempts at seeing quantum effects in micron to nanometre 
scale mechanical systems. Experiments coupling such os- 
cillators to mesoscopic electronic position detectors have 
seen evidence of quantum back-action and back-action 
cooling [T], and have demonstrated continuous position 
detection at a level near the fundamental limit placed 
by quantum backaction [2j [3l H] . Complementary to this 
work, experiments using optomechanical systems (e.g. a 
cantilever coupled to a optical cavity) have been able 
to cool micromechanical resonators by several orders of 
magnitude, using either passive O [6l |7l [8] or active 
(i.e. feedback-based) approaches [3] . 

Despite these recent successes, seeing truly quantum 
behaviour in a mechanical resonator remains a difficult 
challenge. If one is only doing linear position detection, 
the quantum behaviour of an oscillator is almost per- 
fectly masked. Non- linear detector-oscillator couplings 
allow one to probe quantum behaviour such as energy 
quantization [101 IHl IH GSl E]; however, generating 
such couplings is generally not an easy task. Quan- 
tum behaviour could also be revealed by coupling the 
resonator to a qubit [TSl [111 IZ]; this too is challeng- 
ing, as it requires relatively large couplings and a highly 
phase coherent qubit. Here, we consider an alternate 
route to seeing quantum behaviour in a mechanical oscil- 
lator, one that requires no qubit and only a linear cou- 
pling to position. As was first suggested by Braginsky 
and co-workers [18^ il9i , by using an appropriately driven 
electromagnetic cavity which is parametrically coupled 
to a cantilever, one can make a measurement of just a 
single quadrature of the cantilever's motion. As a result, 
quantum mechanical back-action need not place a limit 
on the measurement accuracy, as the back-action affects 



only the unmeasured quadrature. One can then make 
(in principle) a perfect measurement of one quadrature 
of the oscillator's motion. This is in itself useful, as it al- 
lows for the possibility of ultra-sensitive force detection 
[5ni[2I]. Perhaps even more interesting, one expects that 
such a measurement can result in a quantum squeezed 
state of the oscillator, where the uncertainty of the mea- 
sured quadrature drops below its zero point value. 

While the original proposal by Braginsky is quite old, 
there nonetheless does not exist a fully quantum theory of 
the noise and back-action of this scheme; moreover, there 
exists no treatment of the measurement- induced squeez- 
ing. In this paper, we remedy this situation, and present 
a fully quantum theory of measurement in this system. 
We calculate the full noise in the homodyned output sig- 
nal from the cavity (an experimentally measurable quan- 
tity), and derive simple but precise conditions that are 
needed to beat the conventional quantum limit on the 
added noise of a position detector [THl [201 [HI [H] ■ Us- 
ing a conditional measurement approach, we also discuss 
the conditions required to squeeze the mechanical res- 
onator, and demonstrate how feedback may be used to 
unambiguously detect this squeezing. Our results are es- 
pecially timely, given the recent experimental successes in 
realizing cavity position detectors using both supercon- 
ducting stripline resonators [SJ as well as optical cavi- 
ties [3 [6l [71 [3 [25] ; our theory is apphcable to both these 
classes of systems. Note that Ruskov ct al. [2ff recently 
analyzed a somewhat related scheme involving strobo- 
scopic measurement of an oscillator with a quantum point 
contact. Unlike that scheme, the system analyzed here 
should be much easier to implement, being directly re- 
lated to existing experimental setups; our scheme also 
has the benefit of allowing significant squeezing without 
the need to generate extremely fast pulses. 

The remainder of this paper is organized as follows. In 
Sec. [IT] we give a heuristic description of how one may 
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realize back-action free single-quadrature detection, in- 
troduce the Braginsky two-sideband scheme, and give a 
synopsis of our main findings. In Se c. |III| we provide the 
details of our calculations, and Sec. IIVI concludes. 



II. MODEL AND MAIN RESULTS 

A. Basic idea behind single quadrature detection 

Consider a high-Q mechanical oscillator having fre- 
quency ojm and annihilation operator c. We will use X 
and Y to denote the cosine and sine quadratures of the 
oscillator's motion. Using Schrodinger operators c and 
c\ the operators associated with the quadratures are: 



X 



1 

71 



Y = 



(la) 
(lb) 



The Heisenberg-picture position operator x(t) is then 
given by the Heisenberg picture operators X{t) and Y{t) 
via: 

x{t) = v^Xzpt (^X{t) ■ cos wm^ + Y^t) sinujMt^ (2) 

as expected. Note that X and Y are canonically conju- 
gate: 



X.Y 



(3) 



Also note that the definition of the quadrature operators 
relies on having an external clock in the system which 
defines the zero of time. 

In general, X(t) and Y(t) will vary slowly in time (in 
comparison to lum) due to the external forces acting on 
the oscillator. Our goal will be to make a weak, con- 
tinuous measurement of only X(t), using the usual kind 
of setup where the position of the oscillator is linearly 
coupled to a detector. We will use a detector - oscillator 
coupling Hamiltonian of the form: 



ifint = -AxF 



(4) 



where F is some detector operator. It represents the 
force exerted by the detector on the oscillator; in the 
cavity-position detector we will consider, F will be the 
number of photons in an electromagnetic cavity. Ideally, 
if we only measure X, the back-action of the measure- 
ment will only affect the unmeasured quadrature Y, and 
will not affect the evolution of X at later times. The hope 
thus exists of being able to make a back-action free mea- 
surement, one which is not subject to the usual standard 
quantum limit [20l [22l |23] . 

As is discussed extensively in Refs. [20| and |2T] single 
quadrature detection with the interaction Hamiltonian in 
Eq. can be accomplished by simply modulating the 
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FIG. 1: Schematic picture of the setup studied in the text. 
A cavity is driven by an input beam that is amplitude- 
modulated at the mechanical frequency ujm of the movable 
end-mirror. The radiation pressure force, as well as the cav- 
ity and mechanical frequencies and decay rates are indicated. 



coupling strength A at the oscillator frequency. Setting 
A = A{t) = 2{A/ x^pt) cos{ujMt) , Hint becomes: 



int 



= -V2AF- 

X{1 + cos{2uJMt)) + Ysm{2ujMt) 



(5) 



In a time-averaged sense, we see the detector is only cou- 
pled to the X quadrature; we thus might expect that the 
(time-averaged) output of the detector will tell us only 
about X. In principle, this in itself does not imply a lack 
of back-action: via the coupling to Y, noise in F could af- 
fect the dynamics of X. To prevent this, we need the fur- 
ther requirement that the detector force has no frequency 
components near ±2u;m- In this case, the effective back- 
action force F siii(2LUMt) will have no Fourier weight in 
the narrow-bandwidth around zero frequency to which X 
is sensitive, and it will not affect X. Note that Ruskov 
et al. [26j recently considered a linear position detection 
scheme where the effective coupling constant is harmon- 
ically modulated; however, their scheme does not satisfy 
the second requirement above of having a narrow-band 
back-action force. 



B. Model 

We now consider a specific and experimentally- 
realizable system which can realize the above ideas; this 
system was first proposed by Braginsky [181 US- As 
shown schematically in Fig. [l] the setup consists of a 
high-Q mechanical oscillator which is parametrically cou- 
pled with strength ^ to a driven electromagnetic cavity: 



H 



h{uji^ — Ax) [a^ a — (a^a)] 



-H 



M ■ 



(6) 



where ujr is the cavity resonance frequency, Hm = 
hu)M<^ c is the mechanical oscillator Hamiltonian, HdrWc 
describes the cavity drive, H^. describes the cavity damp- 
ing, and is the mechanical damping. Note this same 
system was recently shown (in a similar parameter regime 
to what we will require) to allow back-action cooling to 
the ground state [57]. 
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Assuming a one-sided cavity, standard input-output 
theory [28l|29] yields the Heisenberg equation of motion: 



a = I —ilor 



a - ^/Kbin{t). 



(7) 



where k is the cavity damping, and bin describes both the 
drive apphed to the cavity as well as the noise (quantum 
and thermal) entering the vacuum port. 

To implement back-action evasion in the cavity system, 
we will consider the case where w/j ^ ujmj and take the 
resolved-sideband or "good cavity" limit, where ujm ^ k. 
We will also take an amplitude-modulated cavity drive of 
the form: 



(Si„(t)) = Kit) 



2 



sin(c^Mi)e-*"«* 



(8) 



The same resolved sideband limit is required to achieve 
ground state cooling [27|; all that is different from the 
setup here is the nature of the drive. Here one drives the 
cavity equally at both sidebands associated with the os- 
cillator motion, while in the cooling case, one only drives 
the red-detuned sideband. 

To proceed, we may write the cavity annihilation oper- 
ator a as the sum of a classical piece a(t) and a quantum 
piece d: 



a{t) = a{t) + d{t) 



(9) 



a{t) is determined solely by the the response of the cavity 
to the (classical) external drive bin{t)- In the long-time 
Hmit, Eq. ^ yields: 

a{t) = a„,axC0s(wA/i-H(5)e-""«* (10) 



with 



Omax — &LO 



4^M 



5 — arctan (k/wm) 



(11a) 
(lib) 



The phase 5 plays no role except to set the definitions 
of the two quadratures X and Y; thus, without loss of 
generality, we will set it to zero. We will also be interested 
in a drive large enough that Omax 3> 1. 

In contrast to a, d, the quantum part of the cavity an- 
nihilation operator, is influenced by both the mechanical 
oscillator and quantum noise associated with the cavity 
dissipation. Making use of the solution for d and the con- 
ditions ojR ^ ^ and keeping only terms which are 
at least order a, the term in the total system Hamiltonian 
coupling the oscillator to the cavity takes an analogous 
form to Eq. ([5| with: 

A = \ (Axzpt) ftmax (12a) 
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(12b) 



Thus, the chosen cavity drive gives us the required 
harmonically- modulated coupling constant: in a time- 
averaged sense, the cavity is only coupled to the X oscil- 
lator quadrature. Further, the second condition outlined 



in Sec. II A is also satisfied: because k <C wm, F has no 
appreciable noise power at frequencies near ±2wm- As 
such, we expect no back-action heating of the X quadra- 
ture in the resolved-sideband limit k/lum 0. We will of 
course consider the effect of a non-zero but small k/lom 
in what follows. Note that electromechanical realizations 
of the system presented here, which use superconducting 
stripline resonators, can easily achieve the resolved side- 
band limit lum ^ K (e.g. the recent experiment of Ref . [Ml 
achieved lum/i^ ~ 5). The resolved sideband limit has 
also recently been reached for an optical mode propagat- 
ing in a microtoroidal cavity that supports mechanical 
vibrations, with lum/k « 20 |25j . 



C. Back-action 

Working in an interaction picture, one can easily de- 
rive Heisenberg equations of motion for the system, and 
solve these in the Fourier domain (c.f. Eqs. (43 1). As 



expected, one finds that in the ideal good cavity limit 
{k/lom 0), the measured X quadrature is completely 
unaffected by the coupling to the cavity, while the un- 
measured Y quadrature experiences an extra back-action 
force due to the cavity. For finite k/ldm, there is some 
small additional back-action heating of the X quadra- 
ture. The noise spectral density of the quadrature fiuc- 
tuations are given by: 



Sx{uj) 



Sy{(^) 



dte''^U{X{t),X{0)} 



7/2 



LO^ + (7/2)2 



[1 + 2 {ricq + ribad)] (13a) 



dte''^U{Y{t),Y{0)} 



7/2 



Lo^ + (7/2) 



- [1 -I- 2 (rtcq -t- UBA + J^bad)] 



where 



exp 



M 



(13b) 



(14) 



is the number of thermal quanta in the oscillator, tiba 
parameterizes the back-action heating of the Y quadra- 
ture as an effective increase in rieq! in the relevant limit 
7 <C K one has: 



8i2 _ 2 {Ax^ptY 



K7 



(f^max) 



(15) 



We have assumed here that the there is no thermal noise 
in the cavity drive: it is shot noise-limited. 

Finally, nbad parameterizes the spurious back-action 
heating of X which occurs when one deviates from the 
good-cavity limit; to leading order in k/ujm, it is simply 
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jiven by: 



nhad 



nBA 
32 



K 
UJm 



(16) 



Note that there is no back-action damping of either 
quadrature (see discussion following Eqs. (43)). 



D. Output Spectrum and Beating the SQL 

We assume that a homodyne measurement is made 
of the light leaving the cavity. Using the solution to the 
Heisenberg equations of motion (c.f. Eqs. (43 1 ) and stan- 



dard input-output theory, one can easily find the noise 
spectral density of the homodyne current I{t). The in- 
formation about X(t) will be contained in a bandwidth 
^ 7 ^ K around zero frequency. Thus, focusing on fre- 
quencies K, we have simply: 



Si{uj) ^ 



Sx{uj) 



32^2 



(17) 



Here, G is a gain coefficient proportional to the homo- 
dyne local oscillator amphtude, and So represents added 
noise in the measurement coming from both the cavity 
drive and in the homodyne detection. If both are shot 
noise limited, we simply have Sq = I- We can refer this 
noise back to the oscillator by simply dividing out the 
factor G^: the result is the measured X quadrature fluc- 
tuations: 



S 



X.i 



s(^) 



G2 



Sxiuj) ■ 



32A2 



(18) 



Now, note that in the good cavity limit the spurious heat- 
ing of X described by ribad vanishes. Thus, in this limit, 
the added noise term (second term in Eq. (18l) can be 
made arbitrarily small by increasing the intensity of the 
cavity drive beam (and hence A), without any result- 
ing back-action heating of the measured X quadrature. 
Thus, in the good-cavity limit, there is no back-action im- 
posed limit on how small we can make the added noise of 
the measurement (referred back to the oscillator) . In con- 
trast, for small but non-zero k/ujm, one needs to worry 
about the small residual back-action described by nbad! 
one can still nonetheless beat the standard quantum limit 
in this case, as we now show. 

To compare against the standard quantum limit, con- 
sider S'x,mcas(0): 



S 



X,mGas 



(0) 



7 
2 

7 



(1 + 2ncq + 271bad) + -ZZJ^^i 
oAJ\. 

(1 + 2ncq -t- 2nadd) 



(19) 



In the last line, we have represented both the residual 
back-action nbad and the added noise of the measurement 
as an effective increase in the number of oscillator quanta 
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FIG. 2: Plot of added noise in the single quadrature measure- 
ment (measured as a number of quanta, riadd, c.f. Eq. (20l) 



versus the strength of the measurement (measured in ter ms o f 
the back-action heating of the Y quadrature, tiba, c.f. Eq. 15 1. 



Different curves correspond to different values of k/um and 
So, the noise associated with the homodyne measurement; 
So = 1 corresponds to a shot-noise limited measurement. The 
standard quantum limit of riadd = 0.5 is shown as a horizontal 
dashed line. 



by an amount nadd- The standard quantum limit (which 
applies when both quadratures are measured) yields the 
condition riadd > 1/2 [20] HH [23]. Here, we find: 



n, 



^add — JT-bad + 



«7 c _ "BA / « \ , 1 



(20) 



Thus, if we are in the ideal good-cavity limit {k/ujm ~^ 0) 
and shot-noise limited, beating the standard quantum 
limit on Tiadd requires a coupling strong enough that 
fiBA > 1/8: the Y quadrature fluctuations must be 
heated up by at least an eighth of an oscillator quantum. 

In the more general case where k/u>m is finite, one 
cannot increase the coupling indefinitely, as there is back- 
action on X. One finds that for an optimized coupling 
of: 



the minimum added noise at resonance is given by: 



J^add 



8cj 



M 



(21) 



(22) 



Thus, even for moderately small k,/lom, one can make 
TT-add Smaller than the standard quantum limit value (see 
Fig.|2|. 
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E. Conditional Squeezing 



are 



Given that the double-sideband scheme described here 
can allow for a near perfect measurement of the os- 
cillator X quadrature, one would expect it could lead 
to a squeezed oscillator state, where the uncertainty in 
X drops below the zero point value of 1/2. However, 
Eq. ( 13a I indicates that in the good cavity limit, the fluc- 



tuations of X are completely unaffected by the coupling 
to the cavity detector. To resolve this seeming contra- 
diction, one must consider the conditional aspects of the 
measurement: what is the state of the resonator in a 
particular run of the experiment? In any given run of 
the experiment, the oscillator will indeed be squeezed. 
However, the mean value of X will have some non-zero 
value which is correlated with the noise in the output 
signal. Once one averages over many realizations of the 
experiment, this random motion of (X) appears as extra 
noise, and masks the squeezing, resulting in the result 
of Eq. (13a I. We make these statements precise in what 
follows. 

A rigorous description of the conditional evolution of 
the oscillator in the setup considered here can be devel- 
oped in analogy to Ref. [301 . which considered ordinary 
linear position detection using a cavity. For simplicity, 
we focus on the good cavity limit, where k/ujm ^ 0. 
We flrst deflne the parameter fc, a measure of the rate at 
which the measurement extracts information, as: 



V = Tl (47nBA) 



(23) 



where Uba represents as before the back-action heating 
of the Y quadrature, and ?7 = ^ < 1 represents the ef- 
ficiency of the homodyne detection [t] — 1 corresponds 
to being quantum limited). One has k = l/r-^ca-si where 
Tmcas is the minimum time required to resolve a difference 
in {X) equal to the zero point rms value from the output 
of the detector; as we are interested in weak measure- 
ments, we expect k/uM 1- Note that k = Srjk, where 
k is the usual definition of the strength of the measure- 
ment [SI]. The scaled homodyne output signal may then 
be written [311: 



m^Vl{x{t))+at) 



(24) 



where ^(t) is white Gaussian noise. In a given run of the 
experiment, ^{t) will be correlated with the state of the 
oscillator at times later than t. 

In exact analogy to Ref. [30l a simple description of 
the conditional density matrix is possible in the limit 
where k ^ A. In this limit, the oscillator density ma- 
trix is Gaussian, being fully determined by its means 
X — {X) , Y = (Y) and its second moments Vx = 
{{X')),Vy ~ {{Y')) and C = {{{X,Y}/2)). In the 
interaction picture (i.e. rotating frame at the oscillator 
frequency), the equations for the means (the estimates) 



X = -Ix+VlVx^ 
Y = -'^Y+VICS, 
and for the covariances are 



Vx = 
Vy = 
C = 



liVx - Teq) 



-kC^ + k/{Ar,) - 7(W - Tcq) 
-7C - kVxC 



where: 



T — - 



(25a) 
(25b) 



(25c) 
(25d) 
(25e) 



(26) 



We stress that these equations are almost identical to 
the standard equations for conditional linear position 
detection ED], with the important exception that 
terms corresponding to the bare oscillator Hamiltonian 
are missing. In a sense, the scheme presented here effec- 
tively transforms away the oscillator Hamiltonian. 

To find the amount of squeezing in a particular run of 
the experiment, we simply find the stationary variances 
for the oscillator's Gaussian state. We have: 



Vy = 



nBA 



Vx 
C 



^2 (1 + 2neq) (?7nBA) + 1/4 - 1/2 



4 (?7riBA) 



= 



(27a) 

(27b) 
(27c) 



Note first that the result for Vy is in complete agree- 
ment with the unconditional result of Eq. (13b I: the 



measurement back-action heats the Y quadrature by an 
amount corresponding to riBA quanta. In contrast, we 
find that unlike the unconditional result of Eq. ( 13a I , 



the measurement causes Vx to decrease below its zero- 
coupling value: it is a monotonically decreasing function 
of n^A (see Fig. |3]) . This is the expected measurement- 
induced squeezing. Of particular interest is the minimum 
coupling strength needed to reduce Vx to its zero-point 
value: 



"BA 



' ''cq 

■n 



(28) 



In other words, lowering the X quadrature uncertainty 
from a thermal value of (1/2 + rieq) to the ground state 
value of (1/2) requires that we at least increase the Y- 
quadrature uncertainty by the same amount. This mini- 
mum amount is only achieved for a quantum-limited de- 
tector ?7 = 1. 

The equation describing the fiuctuations of the mean 
quadrature amplitudes X,Y can also be easily solved. 
Assuming that X — Y — Q dX the initial time, one al- 
ways has {X{t)) = {Y{t)) = 0, where the average here 
is over many runs of the experiment. In the stationary 
state (i.e. once the variances Vx, Vy and C have attained 
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FIG. 3: Plot of the conditional X-quadrature variance Vx 



(c.f. Eq. (27b|) as a function of the measurement strength 



(parametrized in terms of the back-action heating of the Y 
quadrature, neA, c.f. Eq. (15l); one clearly sees that the X- 
quadrature can be squeezed. Different curves correspond to 
different values of the bath temperature (parameterized by 



c.f. Eq. (14l) and measurement efficiency 7). The solid 



red curve corresponds to n, 



curve to rie 



and Tj = 1\ the dashed black 
f and r] = 1; the dashed-dot blue curve to 



rieq = 0, ?7 = 0.1. The horizontal dashed line corresponds to 
the ground state value of the variance, Vx = 0.5. 



their stationary value), one finds Y{t) ~ with no fiuc- 
tuations. X(t) continues to fluctuate, with an autocor- 
relation function: 



{mm) = 



y|e-7l*l/2 



(29) 



Again the average here is over many runs of the experi- 
ment. 

W e m ay now combine the results of Eq. ( 29 1 , and 
Eq. (25c) to find Vx,tot, the total (unconditionalJX vari- 
ance. One finds the simple result (valid in the stationary 
regime): 



X,tot 



Vx 



Vx 



7^-2 



(30) 



This shows that, as expected, averaging the results of the 
conditional theory over many measurement runs repro- 
duces the result of the unconditional theory (i.e. the fluc- 
tuations of the measured X quadrature are completely 
unaffected by the measurement). 



F. Feedback for true squeezing 

In the previous section we saw how the state of the res- 
onator, once conditioned on the measurement record, is 



squeezed. We can use feedback control to turn this con- 
ditional squeezing into "real" squeezing of the resonator, 
where the full, unconditional oscillator variance Vx.tot 
(c.f Eq. (30)) drops below the zero-point value. This is 



accomplished by applying a time-dependent force to the 
resonator which is proportional to X{t), the measured 
value of the X quadrature. Such a force can be used to 
suppress the fluctuations in the mean value of X, and 
in the limit of strong feedback, can remove them com- 
pletely. The only fluctuations that remain are quantified 
by the conditional variances, which are squeezed. Note 
that a similar approach was considered in Ref. |26j . 
More precisely, if one makes the measurement at rate k 



described above (c.f. Eq. (23)), and applies the feedback 
force F(t) = a'yXsmiijMt in the laboratory frame, the 
result is an effective damping of the X quadrature at a 
rate 07/2. Calculating the fluctuations of the X quadra- 
ture under this feedback (the details of which are given 
in the next section), we find that the total unconditioned 
X quadrature variance reaches a stationary state: 



{l + a+^Vx)Vx 
\ + a 



(31) 



Here, Vx is the conditioned variance given in Eq. (27b I. 
Note that when a — > 0, we again get the result of 
Eq. ( 30 1 : the unconditioned X quadrature variance is 
not affected by the measurement. In contrast, in the 
limit of large a, one has V^^^^ — > Vx- Thus, as claimed 
above, in the limit of strong feedback, the total fluctua- 
tions of the X quadrature are reduced to the conditional 
variance; it may thus be squeezed. 

It is also important to ask how this squeezing will man- 
ifest itself in the measurement signal. Calculating S^^{ljj), 
the spectrum of the homodyne current in the presence of 
feedback, and referring it back to the X quadrature, we 
find: 



S 



1 

k 



(neq + 1/2) + aVx 



(32) 



In the absence of feedback (a — > 0), the output spectrum 
consists of the white added noise of the measurement 
(1/fc) plus the measurement-independent X-quadrature 
fiuctuation spectrum S'x(w); this is in complete agree- 
ment with the unconditional theory (c.f. Eq. ( 19 1) . When 
feedback is turned on, the second Lorentzian term in 
the spectrum is modified; this corresponds to the added 
damping and noise caused by the feedback. In the limit of 
strong feedback, where a ^ 00, we find somewhat sur- 
prisingly that all signatures of the oscillator disappear: 
Sx^raca.si'^) ~^ V^- Thus, in the limit of strong feedback, 
while the fluctuations in the X quadrature are Vx, they 
do not appear at all in the output signal! To understand 
why this is, note that in driving the resonator with a 
force proportional X, we are driving it with a signal that 
is correlated with the noise in the output signal. Thus, 
feedback leads to new correlations between the fluctua- 
tions of X and the output noise. In the limit of strong 
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feedback, one finds that the fluctuations of X have a 
variance Vx, but are perfectly negatively correlated with 
the output noise. The result is that the output noise is 
completely independent of Vx ■ 

The above effect, in which the fluctuations of X vanish 
in the output signal, may be regarded as an example of 
noise squashing [321 1331 [M] , This is when one uses feed- 
back specifically to reduce the fluctuations in the out- 
put signal, rather than to reduce the fluctuations of the 
system being measured. This is possible only because 
the feedback uses the output signal, and thus correlates 
the system's fluctuations with the output fluctuations. 
Strictly speaking, squashing refers only to the output 
signal that is part of the feedback loop (the so-called "in- 
loop signal"), not to the actual system being measured. 

While the existence of squashing in no way invalidates 
the real squeezing produced by the feedback, it does make 
it more difficult to observe this squeezing at the detec- 



tor output. Further, any experimental results may be 
subject to the accusation that the feedback protocol may 
have been incorrectly designed to produce only squash- 
ing, with the result that squeezing of the resonator could 
not be inferred from the spectrum of the output. 

A solution to this problem is to make a second mea- 
surement of the mechanical resonator's X quadrature 
(e.g. by using a second cavity coupled to the resonator). 
The measurement signal from this second measurement, 
l2{t), is not subject to squashing because it is not part of 
the feedback loop. As a result, the measurement noise in 
l2{t) is completely uncorrelated with the feedback signal. 
Since the second measurement is also a QND measure- 
ment of the X quadrature, it does not affect the results 
for Vxtot or 'S'x moas('^) derived above. If the rate of the 
second measurement is A, then the spectrum of its output 
(again, referred back to the oscillator) is 



-ffb 

-^X,mcas 



.2M 



1 

A 



(neq + 1/2) 
{2ojhf + {l + ay 



A 

{2ojhf+(\ + 2~kVxh 



(33) 



where: 



A 



(2neq + !)(! + 2kVxh + «) " c^Vx 



a{2 



2{2n. 



cq 



l)fc 



(34) 



The first term in the spectrum above represents the 
added noise of the measurement (e.g. shot noise), while 
the terms in square brackets are a direct measure of 
the oscillator's X quadrature fluctuations. We see 
that with feedback, these are described by the sum of 
two Lorentzians. The integral of the area under these 
peaks directly yields V^^^^, the total (unconditional) X 
quadrature variance in the presence of feedback: 



2^ 



(35) 



It thus follows from Eq.(31l that in the strong feedback 
limit (a 00), the area under the resonant peak in the 
output spectrum directly yields Vx, and hence a direct 
measure of squeezing. Note that in this limit only the 
second Lorentzian term in S^^^^^2{'^) survives, as one 
has: 



lim A = (2nc 

a — >oo 



+ l)-Vx^Vx(l + 2~kVx/j] (36) 



where we have made use of Eq. (30 1. Thus, for strong 



feedback, the squeezing of the oscillator can now be un- 
ambiguously detected in the output signal of the second 
measurement: one obtains a simple Lorentzian resonance 
whose area is simply Vx- We remind the reader that in 



the same limit, the spectrum of the flrst measurement 
shows no signature of the oscillator. Note that in prac- 
tice, the limit of strong feedback is already achieved when 
a 3> max(l, fc/7 = 4nBA)- 



III. DETAILS OF CALCULATIONS 
A. Spectrum of the Detector Output 

1. Equations of motion 

The Heisenberg equations of motion (in the rotating 
frame) follow directly from Hq and the dissipative terms 
in the total Hamiltonian: 



d - \//^|(i)e'"^* ~iA[c{l + . 



-2iijj Mt\ 



h.c. 



= --d-vn^ty^''' 



-iV2A 

1 - 



X{1 + cos{2ujMt)) +Ysm{2uJMt) (37a) 



3*'^*^* -ii(l-he 



2iujM 



*) (d + d^) 

(37b) 



Here, ^ describes noise in the cavity input operator 
In the limit where there is only quantum noise (i.e. shot 
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noise) in the cavity drive, we have: 



Note from Eqs. (43 1 that there is no back-action damp- 



(38a) 
(38b) 



In contrast, f} describes equilibrium noise due to the 
intrinsic damping of the mechanical oscillator. One has: 



{fi^it)-f)it')) n,^5{t-t') (39a) 
{f,{t) ■ {t')) = {n,^+l)5{t-t') (39b) 

where rioq is a Bose- Einstein occupation factor evaluated 
at energy hujM and temperature Tbath- 

The equations of motion are easily solved by first writ- 
ing them in terms of the quadrature operators X and Y , 
and then Fourier transforming. To present these solu- 
tions, we first introduce the cavity and mechanical oscil- 
lator susceptibilities as: 



xmH = 



1 



1 



(40a) 
(40b) 



and define the back-action force fsA via 
/ba(cj) = ~AV2kxr{u:) (i{uj + ujr) + |(w - ujr)) .(41) 

Note that while ^ describes white noise, the cavity sus- 
ceptibiUty xr{^) ensures that fsAi'-^) is only significant 
around a narrow bandwidth centered about zero fre- 
quency. Note also that we define Fourier transformed 
operators via: 



i(c.) 



— OO 

oo 



t 



dt 



A{t) 



(42a) 
(42b) 



As such, one has 



A{uj) = A\-Lo). 

The solutions of the Fourier-transformed quadrature 
operators then read: 



y'^ (f}(a; -I- ujm) + ifi^ - ^m)) 



+ 



2i 



(43a) 



^(w) = Xm(w) 

-/ba(w) 



i^J^ (fi{uj + ujm) - f]\ijJ ~ ujm)) 

fBA{^ + ZcJm) + fBA{^ - 2Wj\/) 



ing of either quadrature, even when one deviates from 
the good cavity limit by having k/ujm > 0. This is easy 
to understand on a purely classical level. Note first that 
that it is only the cosine quadrature (i.e. d -f d^) of the 
cavity which couples to the mechanical resonator. As the 
cavity is itself a harmonic oscillator, this means that only 
the cavity sine quadrature (i.e. d — c?^) will be affected 
by the resonator motion. As the cavity cosine quadra- 
ture provides the back-action force on the resonator (c.f. 
Eq. (12b)), it thus follows that the back-action force is 



completely independent of both quadratures of the me- 
chanical resonator's motion. There is thus no back-action 
damping, as such damping requires a back-action force 
which responds (with some time-lag) to the motion of the 
oscillator. 

Equations ( 13 1 for the noise spectra of X and Y at 



frequencies to <C n now follow directly from Eqs. (43 1 



and Eqs. (38),(39l which determine the noises ^ and 1). 



2. Output Spectrum and Beating the SQL 

Standard input-output theory [351 122] yields the fol- 
lowing relation between &outj the field leaving the cavity, 
and 6in, the field entering the cavity: 



bout{t) = bin{t) + VKa{t) 



(44) 



In our case of a one-sided cavity, this relation becomes in 
the lab (i.e. non-rotating) frame: 



6out('^) = bout{'^) + 



-i{uj — ujji) — k/2 
-i{ijj — ojr) + k/2 



-'iAV2KXR{i-J ~ uJr) ■ X{uj — ujr) (45) 

The first term on the RHS simply represents the out- 
put field from the cavity in the absence of the mechan- 
ical oscillator and any fluctuations. It will yield sharp 
peaks at the two sidebands associated with the drive, 
uj = ujR± ujm- The second term on the RHS of Eq. (45 1 



represents the reflected noise of the incident cavity drive. 
This noise will play the role of the "intrinsic output noise" 
or "measurement imprecision" of this detector. 

Finally, the last term on the RHS of Eq. ( 45 1 is the 
amplified signal: it is simply the amplified quadrature 
X of the oscillator. We see that the dynamics of X will 
result in a signal of bandwidth ~ 7 centered at the cav- 
ity resonance frequency. This can be detected by making 
a homodyne measurement of the signal leaving the cav- 
ity. Using a local-oscillator amplitude b^oit) — ?i?e~*'^^* 
with B real, and defining the homodyne current as: 

m={blo{t)+blut{t)) {b Lo{t) +boutit)^ (46) 



(43b) one finds that the fiuctuating part of / is given in 
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frequency-space by: 



(47) 



,iuj + k/2 

iuj — k/2 



^ (eVfl + ^) - ?^^(-^fi + ^)) 



The signal associated with the oscillator will be in a band- 
width ~ 7 ^ k: for these frequencies, the above expres- 
sion simplifies to: 



I{lu) = -B 



AA 



(48) 



Using this equation along with Eqs. (43a I, (38 1 and 



( 39 1 , it is straightforward to obtain the result for the 



homodyne spectrum Si{uj) given in Eq. (17 1. 



B. Conditional Evolution 

To derive the stochastic master equation describing the 
conditional evolution of the resonator under the double 
sideband measurement scheme, (that is, the evolution 
given the continuous stream of information obtained by 
the observer), one uses a procedure that is essentially 
identical to that given in Ref. [30] ■ Once we have moved 
into the interaction picture (in which the quadratures are 
QND observables) , the displacement picture |3S] (that is, 
separated d into d and d as per Eq. and made the 
rotating-wave approximation, the Hamiltonian for the 
combined cavity and resonator system is 



H = -V2A{d + d'')X. 



(49) 



We now perform homodyne detection of output from the 
(one-sided) cavity, and as a result the evolution of the 
system is given by the quantum optical stochastic master 
equation [35| l36] 



da= -^[H,a]dt + KV[d\adt + y^n[~id]adW, (50) 

where a is the joint density matrix of the two systems 
as before rj is the detection efficiency, and k is the cavity 
decay rate. The superoperators V and Ti. are given by 

2D[c]cr = 2cacc' - c'ca - ac'c, (51) 
H[c]a ^ ca + ac^ -Ti-[c(T + (jc^](j, (52) 

for an arbitrary operator c. 

We now wish to obtain an equation for the evolution 
of the resonator alone. This is possible so long as the 
cavity decay rate is fast compared to the timescale of the 
cavity-resonator interaction. That is. 



e < 1, 



(53) 



This means that the light ouput from the cavity spends 
sufficiently little time in the cavity that it continually 
provides up-to-the-minute information about the oscilla- 
tor. With this large damping rate, the fluctuations of the 
light in the cavity about the average value d are small, 
and we can thus expand the cavity state described by the 
operator d about the vacuum: 

<y = Poo|0)(0| + (pio|l)(0|+H.c.) 

+pti|l)(l| + (P2o|2)(0| +H.C.) + 0(e3). (54) 

The density matrix for the resonator is then given by 

p = Trc[a] =poo + Pii + 0(e'). (55) 

where Trc denotes the trace over the cavity mode. From 
the master equation (Eq.(50l) we then derive the equa- 
tions of motion for the pij . Adiabatic elmination of the 
off-diagonal elements poi Sind po2 (described in detail 
in Ref. [3ir) allows us to write a closed set of equations 
for the diagonal elements poo and pn. The result is a 
stochastic master equation for p = pOO + pll, which is 

dp = k[X,[X,p]]dt+ y/2iikn[X]pdW, (56) 

where the measurement strength k = AA^ /k. Defining 
k — 8r]k, and making a Gaussian ansatz for the quan- 
tum state, we find Eqs. (25a I - (25el for the means and 



variances of the quadratures X and Y . 



C. Squeezing via Feedback Control 

There are three formulations that can be used to ana- 
lyze the behavior of an observed linear quantum system: 
the Heisenberg picture (the input-ouput formalism), the 
Schrodinger picture (the SME) and the equivalent classi- 
cal formulation, introduced in Ref. [30^. We have already 
used the first two methods in our analysis above. To an- 
alyze the effect of feedback we now use the third. The 
equivalent classical formulation is given by the equations 



dx = -^xdt + \/ jTcqdWx 



(57) 



dy = —ydt + ^-fT^^dWy + ^dVi + ^dVib%) 



\fk 



where x and y are now classical dynamical variables, and 
as always the noise sources, Wi and Vi, are mutual un- 
correlated Wiener processes. We have now included two 
measurements of the x quadrature, one with strength 
k and the other with strength A, for reasons that will 
be explained below. The measurement records (i.e. the 
homodyned output signals) for these measurements are 
given by 



dli = Vhxdt + dUi 
dl2 = \/\xdt + dU2- 



(59) 
(60) 
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Once again the Ui are mutually uncorrelated Wiener pro- 
cesses. Of interest are the quantities and {X)2, 
which are (respectively) the two observers' estimates of 
the X quadrature. Note that these are not the same as 
X above. When A = 0, so that there is no second mea- 
surement, the equation of motion for {X)i is naturally 
that given by Eqs (25a I - (25el. With the second mea- 



surement, the dynamics of the means and variances for 
the first observer become 



dXi = 


-{-//2)Xdt+ VkVxdUi 


(61a) 


dYi = 


-i-f/2)Ydt+ VlcdtJi 


(61b) 


Vx = 


-fcy| - ^{Vx - Tcq) 


(61c) 


Vy = 


-kC^ + 2k + 2X- 7(W - Tcq) 


(61d) 


C = 


-jC - kVxC 


(61e) 



where k — fc/(8ryi) and A = A/ (8772) are the strengths of 
the respective measurements (under the usual definition 
of measurement strength |31|), and the rji are the respec- 
tive efficiencies of the measurements. We also introduce 
a fourth set of noises, Ui, where Ui appears in the above 
equations for the first observer, and U2 would appear in 
the equations for the second observer, although we will 
not need those here. The Ui are given by [3D] 

dU, = d/, - X,dt = Vl{x - Xi)dt + dU^. (62) 

While it is not obvious, it turns out that the dUi are 
also mutually uncorrelated, and uncorrelated with all the 
other noise sources. 

Armed with the above equations, we now introduce 
feedback into the system. We apply a continuous feed- 
back force F{t) = a'^ X sm{LO Mt) to the system in the 
lab frame. Discarding rapidly oscillating terms (making 
a rotating-wave approximation), this results in the fol- 
lowing dynamics for the system: 



dx = 
dy = 



^Xi )dt^ y ,^cq 
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ydt + J jT,qdWy 



+ dk/AdVi + VA/W2 



7feqdVK, (63) 



(64) 



The feedback provides a damping force on the X quadra- 
ture with a rate aj/2. Since applying a known force to 
the system cannot change the observers' uncertainty re- 
garding the classical coordinates, the equations of motion 
for the variances for both observers remain the same. The 
equations of motion for the means however also pick up 
exactly the same damping terms. Thus for observer one 
we have 

li = -^-^^^X, + VlVxUi (65) 
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'kCUi. 



We now want to calculate the variance of the X 
quadrature under this feedback protocol, and also the 
spectrum of the output signal for both observers. Since 
the X and Y quadratures are not coupled, we need 
merely solve the two coupled equations 



7 



^1 - - 



017 

0:7 



Xi + ^iT,^W^ 
+ kVx)Xi + kVxx 



(67) 

VlVxUi, (68) 



where we have used Eq.(62) to write the equation for 



Xi in terms of Ui rather than Ui. The unconditional 
variance of the X quadrature under feedback, which we 
will denote by Vxtot^ given by the variance of x. 
Solving for the steady-state value of V^j^^. using the 
usual techniques of Ito calculu s, a nd using the fact tha t 
foq = Vx + {kh)V'^ (c.f. Eq. (|30|), we obtain Eq. ([si I. 
We see that as the feedback strength a tends to infinity, 
Vx, as claimed above. 
To calculate the spectrum of the output signal for the 
first observer we first transform Eqs ( 67 1 and ( 68 1 to the 



frequency domain and solve them. The solution is of the 
form 



x{uj) 
Xliu;) 



M{uj 




7feqVF,(c^) 

kVxM^) 



with 



{W^{u:)W^{u')) = {Ui{uj)Ui{J)) = 5{uj + uj'). 
The output signal for the first measurement is: 

Ii{uo) = \/lx{uo) + Ui{uj). 
The corresponding output spectrum is defined via: 

{h{Lo)h{oj'))^Si,{uj)5{u + J) 



(69) 



(70) 



(71) 



(72) 



Using Eqs. (69l,(^70|, one finds the zero-frequency spec- 



trum to be given by: 



{l + a + 2kVx/-lf 



(73) 



where we have made use of Eq. (30 1. Referring this back 
to the oscillator by dividing by k results in Eq. ( 32 ) . 

Similarly, the output of the second measurement is 
given by 



hiu}) = V\x{uj) + U2{lo), 



(74) 



Using this definition and Eqs. (69 1, (70 1, and making 



(66) Eq. (ISSl) 



use of Eq. ( 30 ) , we find the output spectrum given in 
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IV. CONCLUSIONS 

In this paper, we have provided a thorough and fully 
quantum treatment of back-action evasion using a driven 
electromagnetic cavity which is parametrically coupled to 
a mechanical oscillator. We have considered both the un- 
conditional and conditional aspects of the measurement. 
In particular, we have derived exactly how strong the 
coupling must be to beat the standard quantum limit, 
and to achieve a conditionally squeezed state. We have 
also shown how feedback can be used to generate true 
squeezing, and how this squeezing can be detected using 
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